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A new effect in graphene in the presence of crossed uniform electric and magnetic fields is pre- 
dicted. Landau levels are shown to be modified in an unexpected fashion by the electric field, leading 
to a collapse of the spectrum, when the value of electric to magnetic field ratio exceeds a certain 
critical value. Our theoretical results, strikingly different from the standard 2d electron gas, are 
explained using a 'Lorentz boost', and as an 'instability of a relativistic quantum field vacuum', ft 
is a remarkable case of emergent relativistic type phenomena in non-relativistic graphene. We also 
discuss few possible experimental consequence. 



Graphene has been in the forefront of nano electronics 
and quantum condensed matter physics in the last cou- 
ple of years. The mechanical rigidity, being able to peel 
off single graphene layerfl] and ability to make electrical 
contacts has made the system very appealing from build- 
ing devices and experimental points of view One 
of the remarkable discoveries in graphene is the anoma- 
lous 2c? quantized Hall effect. A variety of rich physics 
and anomalous phenomena are tied to a remarkable 'rel- 
ativistic like' spectrum that electron and holes possess in 
graphene [4i] ^] . This has made graphene important and 
interesting from several points of view in physics. 

In the present letter, we investigate the effect of a uni- 
form electric field, applied along the graphene sheet, on 
its already anomalous Landau level spectrum. We find 
that, within the low energy approximation near the Fermi 
surface (Dirac points), the problem can be exactly solved. 
We find strikingly new effects of electric field on the Lan- 
dau levels which is different from the Landau levels of 
standard 2d electron gas. 

We find that the Landau spectrum gets scaled, for a 
given ky quantum number, by an electric field depen- 
dent dimensionless parameter(/? = ;7~b)- As the value of 
this parameter is increased, spacings between the Lan- 
dau levels decreases. This Landau level contraction is 
consequence of electric field induced quantum mechan- 
ical mixing of Landau levels. The entire Landau level 
structure collapses at a critical value of this parameter. 
Further, the 'relativistic' character of the spectrum (with 
Fermi velocity replacing the velocity of light), leads to a 
novel interpretation of our result in terms of relativis- 
tic boosts and the mixing of electric and magnetic fields 
in moving frames of reference. We confirm our analyt- 
ical result by solving the full tight binding model for a 
graphene sheet in the presence of magnetic and electric 
fields numerically. The collapse seen in the low energy 
approximation, is indeed accelerated in the actual tight 
binding model. 

The modified wave functions and energy spectrum will 
have implications on the nature of quantum Hall break 
down. We briefly touch upon this issue at the end, and 
point out how it could be different from the standard 



quantum Hall break down. 

Electronic states of graphene are well described by 
the tight binding hamiltonian for the tt electrons of the 
carbon atoms. In graphene, the carbon atoms form a 
triangular lattice with a basis of two geometrically in- 
equivalent atoms placed ^ apart, where a = 2.456 A" 
is the lattice constant. The overlap integral between 
the nearest carbon atoms is i ~ 2.71eV. We denote 
the triangular lattice sites by — iiSi + 1262, where 
ei = X and 62 = — + are the basis vectors. 
Ciro- {f = 1,2 and a =tj i) represent the electron annihi- 
lation operators with sub-lattice index r and spin index 
a at R;. The hamiltonian is then written as. 
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where 63 = — (ei + 62). The electronic dispersion for 
graphene has two points in a Brillouin zone which sepa- 
rates the positive and negative energy eigenstates. These 
so called Dirac points are ifi.2 = ±^(^x-|-y). The 
dispersion relation in the proximity of the Dirac points 
is linearly proportional to |k| . The low energy modes 
around these points are described by slowly varying fields 
iprTjaCRi) defined as, 

.,Vr'2.(R,) (2) 
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Where , are the Pauli matrices. The effective hamil- 
tonian for the low energy modes is the Dirac hamiltonian. 



H = VF 
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^f^o- are two 



where vp = is the Fermi velocity, 

component field operators where ri{= 1, 2) is the valley 
index, corresponds to two Dirac points and cr(=|,4) is 
the spin index. The spectrum can be obtained by solving 
the one particle equation to get the linear dispersion, 
e(k) = ±?iwi?|k|. In presence of an external magnetic 
field perpendicular to the graphene plane the one particle 
hamiltonian, h = vpa.Tl, where IT = p-l-eA, The energy 
eigenvalues are 



sgn(n)i 



-hvp 
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FIG. 1: Energy eigenvalues e„,fcy, for electrons in graphene 
computed from the tight binding model for a hexagonal lattice 
subjected to a magnetic field B = 27.3 Tesla (or Ic — 20 a 
where a is the triangular lattice spacing) for a system size 
of 600a X 600a. The plot shows e„.ky in the units of t as 
function of ky, where ky is the wavevector in the y direction. 
Two sets of horizontal lines are Landau levels corresponding 
to the two valleys and n — Landau level and the edge states 
are degenerate 

n is the Landau level index, ky — j^l is the quantum 
number corresponding to translation symmetry along y- 
axis, both n and I are integers (we choose Landau gauge 

A(r) = xBy) and Ic ~ ^^'^ magnetic length. 

Unlike the case of the non-relativistic electron in a mag- 
netic field, where the spectrum has a linear dependence 
on the magnetic field and the non-negative integer valued 
Landau level index, the graphene Landau levels have a 
square root dependence on both magnetic field and Lan- 
dau level index. The degeneracy of each level is given by 
the number of magnetic flux quanta passing through the 
sample. The eigenfunctions are, 

where <j)n{S,) are the harmonic oscillator eigen- functions 
and ^= j-{x + l^ky). 

We now consider the above system in the presence of 
a constant electric field in the i-direction. The single 
particle hamiltonian is then given by, 

h = VFOL.TI + leEx (6) 

The Lorentz covariant structure of the hamiltonian, 
with vf playing the role of the speed of light, can be 
used to solve it exactly 0] . It is known from special rela- 
tivity, iivpB > |E|, then we can always boost to a frame 
of reference where the electric field vanishes and the mag- 
netic field is reduced. We can then use the solution in 
Eq. ^ and boost back to get the exact spectrum of the 
hamiltonian in Eq.®. Here the boost transformation 



FIG. 2: Energy eigenvalues e^.i-y for electrons computed for 
the tight binding model for parameters given in Fig[T]and an 
external electric field E applied along a:;-axis, given by the 
parameter /3 = ^j^g — 0.1. The electric field gives a linear 
ky dependence to the bulk Landau levels whereas it gives a 
constant shift to the edge states. The part of solid line labelled 
'bulk Landau level' are n = Landau levels and parallel lines 
above and below them are Landau levels corresponding to 
positive and negative n respectively. Set of points parallel 
to ky labelled 'edge states' are surface states localised at the 
zig-zag boundary. 

amounts to doing a transformation on the space-time co- 
ordinate system. To implement the above procedure, it 
is convenient to work with the manifestly covariant time 
dependent Dirac equation, 

ihr{d^,^i^A^)^{x>')^Q (7) 

Where x'^ — vpt, x^ — x, x'^ — y, 7" = a^, = ia^, 
7^ — —ia^, — qItt. A'^ = (j), the scalar potential, 
— Ax, = Ay and \E'(a;^) is a two component spinor. 
We now apply a Lorentz boost in the y-direction (perpen- 
dicular to the electric field), 

/ \ _ / cosh6i sinh6i \ f x° \ 

\S:^ ) ~ \ sinhd cosh6l ) \ x^ ) ^ ' 

and x^ — x^ . The wave function transforms, ^'(a;^) = 
e2"H\i>(2:^). Applying the above transformations and 
choosing tanh^ = — ^ — B, we can rewrite the Dirac 
equation in Eq.(I7]), 

(^7°9o + 7'9i + j^{d2 + ^Vl^i^)^ ^{in = (9) 

In the boosted coordinates, where |/?| < 1 ,it is a prob- 
lem of a Dirac electron in a (reduced) magnetic field, B — 
B\Jl — . The time component of the 3-momentum in 
the boosted frame, e„ ^ = sgn(n)-y/2|n|^^(l — /3^)4 is 
not the physical energy eigenvalue of our problem. We 



have to apply the inverse boost transformation to obtain 
the spectrum and eigenfunctions of our problem, 

e„,fc„ = sgn(n)y2H^(l ~ P^Y - hvpPky (10) 
oc e"«e-t-. ( *"(;')*[;.l7;«'' ) (11) 

^-^^(-^'.--<")^) 

The energy eigenvalues of the standard 2d electron 
gas in crossed magnetic and electric fields are given by 
en,ky ^ {n+ Dfiujc - nky% - f The main dif- 

ference between the two besides the \/n and ^fB depen- 
dence, is that the low lying graphene Landau level spacing 
scales as (1 — /J^)', whereas the spacing is independent 
of the electric field in the non-relativistic case. Compar- 
ing the eigenfunctions with and without the electric field 
pn [5]), we see that the effect of the electric field is to 
(un) squeeze the oscillator states as well as to mix the 
particle and hole wave- functions. Squeezing corresponds 
to the change in Ic and the eigenfunctions in Eq. (jlip can 
be expanded as superposition of states in Eq.®. Thus, 
unlike in the usual semiconductor samples, in graphene 
the electric field causes Landau level mixing. Also notice 
in p2p that the location of the gaussian also shifts as a 
function of the Landau level index n, unlike the standard 
2d electron gas. 

As f3 approaches unity, from Eg. p^ we infer that, to 
keep the gaussian shifts within the linear extent of the 
system requires larger values of fcy, which takes us be- 
yond the long wavelength approximation. Moreover the 
Eq. (fTO|) hands a collapse of the Landau level spectrum at 
/3 = 1. One may wonder if the collapse we have found 
is an artifact of the low energy approximation? Interest- 
ingly, we find that in our full tight binding calculation 
the collapse persists, and infact it occurs at a value of /? 
even smaller than unity. 

We have performed extensive numerical computations 
on the tight binding model for graphene with magnetic 
and electric fields, using lattice sizes ranging from 60 x 60 
to 600 X 600. The magnetic field enters through the 
Peierls substitution, t t ^ ^■'^'^ in Eq IJ). A(r) is 
chosen in such a way that the contribution to the phase 
term comes from hopping along one of the three bonds for 
each carbon atom. This enables us to maintain transla- 
tion symmetry along the 62 axis of the triangular lattice. 
The problem then reduces to the ID Harper equation. 

, „, k2a + nnp 

e^l^m = 2tCOS{ ^ )<P2,ni + r02,ni + l 

.k2a + niip 

e(P2,ni = 2tC0S( )4>l,ni + tcpi^m^i (13) 
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FIG. 3: (a) shows the energy eigenvalues around the Dirac 
point plotted as function of ky for (3 = 0.9. The collapse 
can be clearly seen (b) shows the modulus of the eigenvalues 
Ifn.fcyl for a value of ky = 2.785 computed from the tight 
binding model for system size 600a x 600a, magnetic field 
B — 27.3 Testa or Ic — 20 a, electric field given by parameter 
(3 = 0.0, 0.5, 0.9, plotted as a function of n. 



Here tp is the magnetic fiux passing through each plaque- 
tte, k2 is the wave vector and ni is the ei component of 
triangular lattice coordinate. 

We choose the value of the magnetic field such that 
-/j S> S> a, where L is the linear extent of the system. 
The condition Ic^ a ensures that we stay away from the 
Hofstadter butterfly kind of commensurability effects on 
the spectrum and L ^ Ic ensures that a large number 
of cyclotrons orbits fit in the sample. For our numerics, 
we expressed all energies in units of t and all lengths in 
units of a. 

Fig.ll]) shows the results of our numerical investiga- 
tion for zero and Fig.(l2|) for a finite {(3 = 0.1) electric 
fields. Fig.(IT|) shows the spectrum at low energies and 
the eigenvalues that are constant w.r.t ky are the Lan- 
dau levels. They have -^/n behaviour and are in excellent 
agreement with the analytical result and the eigenstates 
that vary with ky are the chiral edge states responsible 
for the quantum Hall current. In our numerics the lat- 
tice has zig-zag edges at the two ends along ei. It is well 
known that the zig-zag edges result in zero-eigenvalue 
states even in the absence of external magnetic field Q. 
These zero-eigenvalues surface states which are localised 
at the boundary and n = Landau level forms degen- 
erate set of states as shown in Fig.([l|). However, Fig.Q 
shows that this degeneracy gets lifted in presence of an 
electric field. For small electric fields the wavefuntion of 
these edge states continues to be localised near edge of 
the sample. Any difference in eigenvalues of these sur- 
face states is because of the potential seen by them due 
to the externally applied electric field at the two edges of 
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the sample. A characteristic feature of these edge states 
is that they don't vary with the wavevector, whereas the 
Landau levels develop a linear ky dependence with elec- 
tric field. 

Fig.lISb), shows scaling of Landau levels for a given 
ky value. For zero electric field we see an excellent match 
between analytics and numerics. And for the case of finite 
electric fields we see a systematic deviation from exact 
results as we suspected from our exact result. As /? ^ 1, 
the tight binding results shows a faster collapse. Fig.([3K) 
shows the collapse has already occurred at /3 = 0.9, near 
one the Dirac points. 

We show below that one of the consequences of the 
Landau level contraction (jlOp and the n dependent guas- 
sian shift is the possibility of a 'dielectric break- 
down', which is different from the conventional ones. The 
single particle spectrum and states we have obtained thus 
far (for a given E and B) can be used to construct stable 
many-body quantum Hall ground states. However, the 
external electric field not only modifies the single particle 
wave function and spectrum, but can also destabilise the 
ground state through spontaneous creation of particle- 
hole pairs; i.e., by a dielectric breakdown. We present 
a simple formula for dielectric breakdown, without giv- 
ing full details. It has an unusual dependence on the 
length scale over which the potential fluctuates and on 
the Landau level index n. This peculiar feature is absent 
in standard quantum 2d electron systems . Specifically 
we find that for slowly varying electric field fluctuations 
over a length scale £e and for large Landau level index, 
the critical voltage for breakdown is given by: 



T4 « l±Kn . 

e \ \£e, 



(14) 



where k is a constant of the order of unity, which depends 
on the strength of the electric field fluctuations and A„ 
is gap between levels n and n + 1. This means that if 
we have an electric field, non-uniform over a nanoscopic 
scale (£_E ~ Ic), it will cause local breakdown even be- 
fore the critical field is reached. Such situations can be 
created through in plane or out of plane charged impu- 
rities or STM tips, in addition to external electric fields. 
It is interesting that such an anomalous local breakdown 
is Landau level index n dependent, we expect that the 
quantum Hall breakdown should be qualitatively differ- 
ent for n = and n ^ within graphene. 

In the light of new spectroscopic experiments Q, we 
claim that the contraction in Landau level spacing and 
the collapse can be observed at fields attainable in labora- 
tories. The gap between n = and rt = 1 for B ~ 1 Tesla 
is ~ 35 meV, for ~ 3 x 10^ Vm^^, 10% reduction in the 
gap is expected. And the collapse of the Landau levels 
should also be observeable by applying E ^ lO^V^m^^. 
In the context of quantum Hall breakdown, the depen- 
dence of critical voltage on £e as given in Eq. (fT4|) suggests 



that the breakdown phenomena should be different from 
what we observe in standard 2d quantum Hall system. 
Moreover graphene's Landau level index dependence on 
Vc, we expect the breakdown phenomena is going to be 
different for n ^ from that of n = 0. 

It will be interesting to study graphene from the point 
of view of the present paper. As quantum Hall phenom- 
ena are beginning to be seen in p yrolytic graphite 1 10] and 
possibly in carbon eggshells Uj , it will be very interest- 
ing to study electric field effects in these systems as well, 
to confirm our predictions. 

In summary, we have made a theoretical prediction 
of a remarkable phenomena in graphene: Landau level 
contraction and an eventual collapse, induced by crossed 
electric fields. The local dielectric breakdown has a pe- 
culiar length and Landau level index dependence. These 
phenomena, not known in the standard 2d electron gas, 
is a consequence of the relativistic type spectrum of low 
energy electrons and holes in graphene. 
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